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ABSTRACT 

Context. The variation of the dimensionless fundamental physical constant yu = mp/nig - the proton to electron mass ratio - can be 
constrained via observation of Lyman and Werner lines of molecular hydrogen in the spectra of damped Lyman alpha systems (DLAs) 
in the line of sight to distant QSOs. 

Aims. Our intention is to max out the possible precision of quasar absorption spectroscopy with regard to the investigation of the vari- 
ation of the proton-to-electron mass-ratio yu. The demand for precision requires an understanding of the errors involved and effective 
techniques to handle present systematic errors. 

Methods. An analysis based on UVES high resolution data sets of QSO 0347-383 and its DLA is put forward and new approaches to 
some of the steps involved in the data analysis are introduced. We apply corrections for the observed offsets between discrete spectra 
and for the first time we find indications for inter-order distortions. 

Results. Drawing on VLT-UVES observations of QSO 0347-383 in 2009 our analysis yields A/i/jj = (4.3 ± 7.2) x 10"*^ at z-^b. = 3.025. 
Conclusions. Current analyzes tend to underestimate the impact of systematic errors. Based on the scatter of the measured redshifts 
and the corresponding low significance of the redshift-sensitivity correlation we estimate the limit of accuracy of line position mea- 
surements to~220ms"', consisting of roughly 150ms"' due to the uncertainty of the absorption line fit and about 150ms"' allocated 
to systematics related to instrumentation and calibration. 

Key words. Cosmology: observations - quasars: absorption lines - quasars: individual: QSO 0347-383 



1. Introduction 

The Standard Model of particle physics contains several funda- 
mental constants whose values cannot be predicted by theory and 
need to be measured through experiments (Fiitzsch 2009 ). They 
are the mass of the elementary particles and the dimensionless 
coupling constants which are assumed time-invariant although in 
theoretical models which seek to unify the four forces of nature 
they vary naturally on cosmological scales. The fine-structure 
constant a = e^ I {Aneffic) and the proton-to-electron mass ratio, 
jU — niplme are two constants that can be probed in the labo- 
ratory as well as in the Universe by means of observations of 
absorption lines due to intervening systems in the spectra of dis- 
tant QSOs and have been the subject of numerous studies. The 
former is related to the electromagnetic force while the latter is 
sensitive primarily to the quantum chromodynamic scale (see, 
i.e., Flambaum l20041 l. 

A probe of the variation of ji could be obtained by compar- 
ing relative frequencies of the electro-vibro-rotational lines of 
H2 as first applied by Varshalovich and Levshakov ("1993) af- 
ter Thompson (1 19751 1 proposed the general approach to utilize 
molecule transitions for /i-determination. The original paper by 
Thompson ( 1975 ) did not take into account the different sensitiv- 
ities within the molecular bands, which is the key of the modern 
approach. 

The method is based on the fact that the wavelengths of 
vibro-rotational lines of molecules depend on the reduced mass. 



M, of the molecule. For molecular hydrogen M - nip/l so that 
the comparison of an observed vibro-rotational spectrum with 
its present analog will give information on the variation of m^f 
and m^. Comparing electro-vibro-rotational lines with different 
sensitivity coefficients gives a measurement of//. 

The observed wavelength /lobs,/ of any given line in an ab- 
sorption system at the redshift z differs from the local rest-frame 
wavelength A^j of the same line in the laboratory according to 
the relation 






(1) 



where Kj is the sensitivity coefficient of the /th component com- 
puted theoretically for the Lyman and Werner bands of the 
H2 molecule (Varshalovich & Levshakov 119931 Varshalovich 
& Potekhin 1995^ Potekhin et al. [T9981 Meshkov et al. l2007l 
Ubachs et al. 2007 ). 

It is useful to measure variations in velocities with compar- 
ison to the redshift of a given system defined by the redshift 



position of the lines with Kj 
redshift ^/: 



AV,. 



(i = 



1+z 



Ki 






0, then introducing the reduced 



(2) 



* Based on observations taken at ESO Paranal Observatory, Program 
083.A.0733 



The velocity shifts of the lines are linearly proportional to Afi/ju 
which can be measured through a regression analysis in the AV,- 
Kj plane. 

This method was used to obtain a bound on the secular vari- 
ation of the electron-to-proton mass ratio at Afi/fi - (-1.8 ± 
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3.8)xlO"^ fromobservationsof the newly discovered H2 absorp- 
tion systems at Zabs =3.0 towards QSO 0347-383 (Levshakov 
et al. 2002 ). Subsequent measures of the absorption systems of 
QSO 0347-383 and QSO 1232+082 provided a strong indica- 
tion of a variation (2.4 + 0.6) x 10"^, i.e. at 3.5 cr (Reinhold 
et al. I2OO6I Ubachs et al. ISoOTl l. Earlier works (Ivanchik et al. 
12005 ^) also find hints for variation, but are still dominated by 
inaccuracies of the laboratory wavelengths. However, more re- 
cently King et al. (2008), Wendt & Rei mers ([2008] ), Tho mpson 
et al. ([2009h ). Wendt and Molaro ( 120111) . King et al. ("20ri ) with 
An/fi = (0.3 + 3.7) X 1 0-^ at Zabs = 2.81 1 towai-ds PKS 0528-250 
and Bagdonaite et al. (2012) with A/z/yK = (-6.8 ±27.8) x lO"*" at 
Zahs - 2.426 towards QSO 2348-011 reported a result in agree- 
ment with no variation. 

The more stringent limits on Afi/fj. have been found from the 
combination of three H2 systems at Afi/fi - (2.6 + 3.0stat)x 10"''( 
King et al. 2008) and taking into account additional transitions 
from deuterated molecular hydrogen (HD) in King et al. (120111 ). 
A fourth system has provided Afi/fi 
10 ''(Malecetal.lMOl). 



(+5.6 ± 5.5stat + 2.9™s) X 



An independent method relies on the inverse spectrum of 
ammonia as shown by Flambaum & Kozlov (2007). Ammonia 
NH3 inversion transitions are very sensitive to changes in fi due 
to a tunneling effect. The sensitivity coefficient of the inversion 
transition can be almost two orders of magnitude more sensitive 
to /i- variation than H2 molecular rotational frequencies. Thus by 
comparing the inversion frequency of NH3(1,1) with a rotational 
frequency of another co-spatial molecule it is possible to bind a 
variation of yu. 



Flambaum & Kozlov ( 120071 ) combine three detected NH3 ab- 
sorption spectra from B0218+357 with rotational spectra of CO, 
HCO+, and HCN to place a Hmit of (0.6± 1 .9) x 10"^ for a look- 
back time of 6 Gyr (redshift z = 0.68). Murphy et al. (2008) with 
newly obtained high signal-to-noise rotational spectra of HCO^ 
and HCN obtained < 1.8 x 10"^ at a 95% CL. 

Henkel et al. (2009) obtained a firm upper limit of 10"^ for a 
look-back time of 7 Gyr (z=0.89) towards PKS 1830-211. This 
method is limited to low redshifts due to the small number of 
NH3 sources in general and to the large line widths and chem- 
ical segregation of different molecules at higher redshifts. For 
sources in the local Milky Way, however, a very strict upper 
limit of \Afj./fi\ < 3 x 10~^ was found utilizing both the ammonia 
method (Levshakov et al. I2010al and Levshakov et al. 1201 Obi ), 
and the methanol method (Levshakov et al. 120 111) . 

Laboratory experiments by comparing the rates between 
clocks based on hyperfine transitions in atoms with a different 
dependence on /i restrict the time-dependence of ft at the level of 
(fi/n),, = (1.6 ± 1.7) X 10-'5 yr-i (Blatt et al. 2008 ). 

In the following we will concentrate on the H2 system ob- 
served towards QSO 0347-383 to trace the proton-to-electron 
mass ration fi at high redshift (z^bs - 3.025). The motivation for 
re-analysis of QSO 0347-383 is given by the dramatically en- 
hanced quality of the recent data of this quasar for the purpose 
of setting constraints on Afi/fj.. The single velocity component in 
H2 absorption renders QSO 347-383 a prime target to further in- 
vestigate the impact of wavelength calibration issues. Trying to 
reach a sensitivity of few parts per million everything becomes 
important and the special requirements of the observations as 
described in the following section are absolutely mandatory. 



2. Data 

2.1. Observations 

The recent observations of QSO 0347-383 were performed with 
UVES on VLT on the nights of September 20-24 2009. The jour- 
nal of these observations is given in Table [T] The DlC 2 setting 
was used with blue setting and the 437 nm grating. The CCDs 
were not binned with pixel size of 0.013 - 0.015 A, or 1.12 
km s"' at 400 nm along the dispersion direction. The observa- 
tions are comprised of 1 1 exposures on four successive nights, 
of which 10 exposures were of 5400 s and one of 3812 s. Eight 
of the spectra with setting 437+760 and three the 437+860 set- 
ting, providing a coverage in the blue spectral ranges between 
373-500 nm. QSO 0347-383 has no flux below 370 nm due to 
the Lyman discontinuity of the Zabs = 3.023 absorption system. 
The slit width was set to 0.7" for all observations providing a 



Table 1. Journal of the observations (2009 data). Before and af- 
ter each spectrum, a 30 s calibration frame was recorded 



No. 


Date 


Time 


A 


Exp[s] 


DIMM[arcsec] 


1 


2009-09-20 


05:05:46 


437 


5400 


1.21(0.29) 


2 


2009-09-20 


08:28:48 


437 


3812 


1.73(0.23) 


3 


2009-09-21 


04:45:51 


437 


5400 


1.15(0.12) 


4 


2009-09-21 


06:18:45 


437 


5400 


1.21(0.18) 


5 


2009-09-21 


07:59:24 


437 


5400 


1.08(0.09) 


6 


2009-09-22 


04:41:37 


437 


5400 


0.97(0.17) 


7 


2009-09-22 


06:14:25 


437 


5400 


1.03(0.13) 


8 


2009-09-22 


07:59:19 


437 


5400 


1.00(0.14) 


9 


2009-09-23 


04:24:05 


437 


5400 


1.36(0.26) 


10 


2009-09-23 


05:56:49 


437 


5400 


1.18(0.23) 


11 


2009-09-23 


07:29:35 


437 


5400 


0.95(0.23) 



mean resolving power of A/AA k 66000. Within each order the 
resolving power varies by about 15-20 % being higher at the 
starting wavelength of each order. The average seeing along the 
exposures as recorded by the DIMM at Paranal is given in the 
last column of Table [T] We note, however, that the actual seeing 
at the UT2 of VLT was significantly better than that recorded by 
the DIMM. 

The 1 1 different spectra and the corresponding co-added data 
(bottom) are shown in a region around L4R1 in FigurefTl 

2.2. Reduction 

The last version of the UVES pipelinaH has been used for data 
reduction. The pipeline first uses a set of five bias to make a 
mean bias free of cosmic ray hits which is then subtracted to 
all the two dimensional format images. A pinhole lamp image 
is used for identifying the precise location of the Echelle orders 
which are curved and tilted upwards. 

The pixel-wavelength conversion of the wavelength calibra- 
tion is done by using the corresponding calibration spectrum. 
Murphy et al. (2008) and Thompson et al. (2009) independently 
showed that the standard Th/Ar line list used in the old UVES 
pipeline was a primary limiting factor The laboratory wave- 
lengths of the calibration spectrum were only given to three dec- 
imal digits (in units of A) and, in many cases, the wavelengths 
were truncated rather than rounded from four decimal places 
(see Murphy et al. 20071 ). 



Version 4.9.5 
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Fig. 1. The 1 1 single spectra and the corresponding co-added data (bot- 
tom) are plotted around the region of L4R1 and L4P1 (vertical lines). 
The stronger H2 features can be merely distinguished in single spectra. 



Thompson et al. ( 12009 l l re-calibrated the wavelength solu- 
tions using the calibration line spectra taken during the observa- 
tions of the QSOs and argued that the new wavelength calibra- 
tion was a key element in their null result. 

The new data UVES pipeline has solved these problems. 

The UVES blue frame comprise 32 orders, from absolute 
number 96 to 124, covering the wavelength range 374-497 nm 
while the molecular lines are spread over 18 UVES echelle or- 
ders (from 106-122) covering the wavelength range 380-440 nm. 

More than 55% of the x 400 ThAr lines in the region were 
used to calibrate the lamp exposure by means of a polynomial 
of the 5th order. Typical residuals of the wavelength calibrations 
were of ~ 0.34 mAor ~ 24 m s"' at 400 nm and were found 
symmetrically distributed around the final wavelength solution 
at all wavelengths. By comparison in Malec et al. ( 120101 ) the 
wavelength calibration residuals have been RMS ~ 80 m s"'. 

In our set of observations, calibration spectra were taken be- 
fore and after the object spectra for each night. Calibration lamps 
taken immediately before and after object observation provides 
accurate monitoring of physical variations. Moreover, the cali- 
bration frames were taken in special mode to avoid automatic 
spectrograph resetting at the start of every exposure. Since Dec 
2001 UVES has implemented an automatic resetting of the Cross 
Disperser encoder positions at the start of each exposure. This 
implementation has been done to have the possibility to use day- 
time ThAr calibration frames for saving night time. If this is ex- 
cellent for standard observations, it is a problem for the mea- 
surement of fundamental constants which requires the best pos- 
sible wavelength calibration. Only calibration spectra that are 
attached as template to the OB allow to take the calibration ex- 
posure in exactly the same instrument setting as the science ex- 
posure. These calibrations can be taken upon user's request in 
addition to the ones from the calibration plan. Thermal -pressure 
changes move in the cross dispersers in different ways, thus in- 



troducing relative shifts between the different spectral ranges in 
different exposures. 

It should be emphasized that this effect has not been taken 
into account in the analysis performed so far on UVES data for 
fi variability. 

There are no measurable temperature changes for the short 
exposures of the calibration lamps but during the much longer 
science exposures the temperature drifts generally by 0.1 K, and 
in two cases the drift is of 0.2 K while in other two there is no 
measurable change. Pressure values are surveyed at the begin- 
ning and end of the exposures and changes range from 0.2 to 
0.8 mbar. The estimates for UVES are of 50 m s^' for AT = 0.3 
K or a AP = 1 mbar (Kaufer et al. .2004J . thus assuring a radial 
velocity stability within ~ 50 m s^ 

Individual spectra are corrected for the motion of the obser- 
vatory about the barycenter of the Earth-Sun system and then 
reduced to vacuum. The velocity component along the direction 
to the object of the barycentric velocity of the observatory was 
calculated using the date and time of the midpoint of the inte- 
gration to minimize the influence of changes. The changes in ra- 
dial velocity during exposure induce a symmetric modification 
of the line profile. The absorption profile is not strictly Gaussian 
(or Voight) anymore but rather slightly squared-shaped (since 
the FWHM of the line is ^ 5 km s 'and the smearing of the line 
by Earth motions of + 40 m s^'the effect is negligible. The line 
shapes remain symmetric in any case and possible changes of 
radial velocities during exposure effects only the quality of the 
fit, it does not influence the measured centroid of an absorption 
line. The wavelength scale was then corrected for this motion so 
that the final wavelengths are vacuum wavelengths as observed 
in a reference frame at rest relative to the barycenter 

The air wavelengths have been transformed into vacuum 
by means of the dispersion formula by Edlen (1966V Drifts in 
the refractive index of air inside the spectrograph between the 
ThAr and quasar exposures will therefore cause miscalibrations. 
According to the Edlen formula for the refractive index of air, 
temperature and atmospheric pressure changes of 1 K and 1 
mbar would cause differential velocity shifts between 370 nm 
and 440 nm of ~ 10 ms"'. 



3. Preprocessing 

3.1. Spectral radial velocity shifts 

Wendt and Molaro (201 1) showed the presence of overall shifts 
between spectra obtained with slit spectrographs such as UVES. 
For checking such a possibility in the new data set we obtained 
a median velocity per each exposure by fitting as many lines as 
possible in each of the eleven spectra. The 50 H2 lines were ini- 
tially selected for that, but due to the relatively low SN spectra 
in single exposures, not all could be fitted, and for instance most 
of the weak J=0 transitions have been missed. The median radial 
velocity in respect to the chosen absorber redshift derived from 
all detected lines within an individual exposure ranges from - 
200 m s^'to + 354 m s^'for the 1 1 spectra. Due to the low qual- 
ity of the fits to single lines of individual spectra with low signal- 
to-noise, the median was chosen. Despite the low quality of the 
individual fits the median velocities listed in Table|2]are well de- 
fined as shown in Figure l2] Bootstrapping is a useful practice of 
estimating properties of an estimator (i.e. its error). For exam- 
ple 48 line positions were determined in spectrum 8 (see Table 
01. Their median velocity offset corresponds to 310 m s '. The 
bootstrap histograms were implemented by constructing a num- 
ber of resamples of the observed line positions (and of equal size 
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Table 2. Median velocities of single spectra 



spectrum median velocity [km s ] lines considered 



1 


0.075 


46 


2 


0.058 


42 


3 


-0.136 


41 


4 


0.223 


43 


5 


0.354 


47 


6 


-0.206 


47 


7 


-0.015 


44 


8 


0.310 


48 


9 


-0.051 


40 


10 


-0.107 


47 


11 


-0.004 


48 



0.046 



to the observed data set), each of which is obtained by random 
sampling with replacement from the original data set. 

For the data at hand the obtained velocity offsets bear no sta- 
tistical significance because of the large scatter of the individual 
position measurements per single spectrum. The resulting offsets 
are, however, comparably well defined which made us confident 
to apply this procedure. The described method had no signifi- 
cant impact on the final result at the current level but it might 
provide a helpful tool in the future to check for potential inter- 
spectra shifts. The resulting velocities reflect the offsets of the 
individual spectra to a reference redshift of the H2 absorber. Of 
the 1 1 shifts a mean velocity offset (relative to the assumed Zabs) 
of 46 m s 'was computecj^ This low residual offset verifies the 
assumed redshift for the absorption system. 




0.1 0.2 0.3 

velocity offset [kms~ 



0.6 



Fig. 2. Exemplary bootstrap histogram of the median position of lines 
in spectrum 8 with respect to z = 3.02489817. The Gaussian fit corre- 
sponds to 302 ± 4 m s"' (compare Tablel2l. 



4. Selection of lines and line fitting 

For the analysis a total set of 50 H2 lines are fitted. This preselec- 
tion is based on earlier line identifications (see Wendt & Molaro 
[2011 ) including curve of growth analysis to avoid blends with 
the Lyman-c forest or other H2 lines. The separate spectra are 



not coadded in this step since the fitting algorithm works on the 
different data sets simultaneously. This allows us to omit further 
rebinning of the data. The fitting code is based on an evolution- 
ary algorithm, which tracks the global minimum via an interac- 
tive process of covering the parameter space (see Quast et al. 
12005b . Each set of fitting parameters is evaluated V\2iX^ in every 
single spectra. Thus, there is no need to redistribute the flux of 
each exposure to pixels of equal central wavelengths. Constant 
velocity offsets, however, potentially influence the shapes of the 
absorption features as would be the case of coadded spectra. 
Evolutionary fitting algorithms are less prone to converge in a 
local minimum rather than find the global minimum. An advan- 
tage over Monte Carlo chain methods as applied for the purpose 
of line fitting in simple cases for example in King et al. (I2008I I is 
the drastically reduced need for computer power. Additionally, 
the principle of evaluating multiple groups of parameters inde- 
pendently allows for consequent parallel computing. 

For each set of lines sharing the same rotational level a com- 
mon column density and a common broadening parameter is 
fitted with simplified pseudo-Voigt-function profiles. For weak 
lines a mere Gaussian profile would suffice since natural line 
broadening has no noticeable impact on the line shapes. In case 
of QSO 0347-383, only a single component is observed in H2 . 
The only free parameter per each individual line is the radial ve- 
locity with respect to the absorber redshift. Out of the 50 H2 lines 
analyzed, eight were excluded since they showed a comparably 
large positioning error of more than 300 m s"'. This is mostly 
due to a continuum highly contaminated by the presence of hy- 
drogen absorption in the environment of the affected lines. 

We note that our procedure is different from the one followed 
by King et al. (I2008I I. They fitted numerous additional compo- 
nents in a region of H2 absorption to narrow down the x^ of the 
fit to the data. 

Figure[3]compares a portion of spectrum used in the analyses 
by Ivanchik et al. (2005), King et al. (I2008I I. Thompson et al. 
(|2009J . and Wendt & Molaro (^CjTTT l with the same portion of 
the new data we are analyzing here. The solid red vertical lines 
mark the H2 component L4R1 and L4P1 and the dotted red lines 
indicate the 12 additional lines in that regiorr] The upper plot 
corresponds to the data of 2009 and reveals that some of the 
extra components clearly recreate the flux observed in 2002 but 
evidently do not correspond to factual properties of the absorber. 

In King et al. (2008) the evolution of ;(f^ with an increasing 
number of additional free lines is the main criterion to fix the to- 
tal number of components. While that approach clearly reduces 
the residuals of the fit, it may not reflect the physical properties 
of the absorber. Though it is likely that the absorber structure is 
too complex to be represented by a single component, we prefer 
to integrate the uncertainty of the true nature of the velocity com- 
ponents into the fitting uncertainty rather than 'generating' com- 
ponents to fill up the flux in a poorly known continuum. Higher 
resolution spectra may verify or falsify some of the decisions on 
additional components and help to distinguish between apparent 
precision (lower ;(f^) and reached accuracy (better description of 
the physical conditions of the absorber) or rather the limit on 
information on the absorber. 



5. Results 

In Figure HI the measured radial velocities of the 42 H2 lines are 
plotted against the sensitivity towards k coefficients of the cor- 



^ A non-zero mean radial velocity directly reflects a deviation from 
the assumed absorber redshift. 



' Their individual positions are extracted from the plot in King et al. 
( I2008l l. 
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Fig. 3. Comparison between the 2009 data {top graph) and, with an 
offset, the original single observation run data of 2002 (9 frames, bot- 
tom) as used by Ivanchik et al. ( 2005 ), King et al. (2008 1, Thompson 
et al. ( 12009b . and Wendt & Molaro (2011 1. The vertical hues indicate 
the positions of the H2 component (solid) and the 12 additional lines fit- 
ted in King et al. ( I2008t . In green the single line fits (L4R1 and L4P1) 
with their corresponding local continuum polynomial fit from this work. 
Note, that the fits for the analysis were not carried out on the rebinned 
and coadded data as shown in this figure. 



responding transition. Any correlation therein would indicate a 
variation of jj at z=3.025 with respect to laboratory values. Table 
[3] lists the broadening parameter and the column density for all 
the lines from one particular J level fit together and consistently 
with one N and one b per rotation level. 



Figure |5] shows the same line data as is Figure |4] but lines 
with similar k values are binned. For a better overview, errorbars 
are omitted. The red data points (crosses) reflect the radial veloc- 
ity within a small sensitivity range (given as x-errorbar). The y- 
errorbars correspond to the standard deviation of the mean value 
of the scatter within such a bin. The scatter within the bins can 
not be attributed to possible variations of fi since it is present 
for basically the same sensitivity parameter. The scatter is of 
the order of 220 m s 'and thus larger than the positioning er- 
ror of the individual lines. That is also reflected by a reduced 
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Fig. 5. The observed 42 H2 lines with their corresponding radial ve- 
locities. In blue, bins of similar sensitivity are plotted. Blue y-errorbars 
reflect the standard deviation of radial velocities within an interval (x- 
errorbar). 
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Fig. 4. Measured radial velocity vs. sensitivity for all 42 lines. Any cor- 
relation would indicate a change in /j. The different symbols and colors 
used correspond to the observed three rotational levels. The errorbars 
reflect the 1 cr errors. 



X of 2.7 for a weighted linear fit to the data (corresponding to 
Afi/H = (1.8 ± 8.2) X 10"'' at za,, = 3.025). The true scatter of the 
data is of the order of 220 m s" ' and constitutes an absolute limit 
of precision. The above mentioned errors of the fitting proce- 
dure require an additional systematic component to explain the 
observed scatter: 



0"obs 



4 



^pos "'" '-^sys' 



with cTobs - 220 m s ^,0-. 



pos 



152 m s 



A direct linear fit to the unweighted data yields: 
An/iu = (4.2 ± 7.7) X 10"^ 



(3) 



and (Tsys = 151 m s '. 



(4) 



Bootstrap analysis is a robust approach to obtain a linear fit 
to the data in Figure [4] and estimate an error based on the true 
scatter of the data. The corresponding bootstrap histogram in- 
cluding Gaussian fit is illustrated in Figure [6] The Gaussian fit 

gives: 



Afi/l^ = (4.3 ± 7.2) X 10" 



(5) 



The data give no hint towards variation of the proton-to- 
electron mass ratio in the course of cosmic time. The uncertain- 
ties in the line positions of the H2 features due to the photon 
noise are estimated by the fitting algorithm. These are shown in 
the errorbars in Figure|4]and reported in Table|4] The mean eiTor 
in the line positioning is of 152 m s"'. Even at first glance the 
given errorbars in Figured appear to be too small to explain the 
observed scatter. 



6. Systematics 

Any process which influences the measured redshift in depen- 
dence with the excitation energy would mimic a variation in jj 
since the transitions from higher excited states naturally show 
a stronger sensitivity towards changes in // (see Varshalovish & 
Levshakov 11993b . Figure IT] marks the wavelength ranges cov- 
ered by the different orders of the CCD spectrum. Decreasing 
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Fig. 6. Bootstrap histogram of 50.000 samples based on all 42 lines. 
The dashed line shows the bootstrap fit of the 2002 data of QSO 0347- 
383. 



wavelengths are plotted rightwards since K sensitivity factors 
are increasing almost linearly with decreasing wavelengths and 
therefore FigurelTlis comparable with Figureffl This new figure 
shows no trend but again a rather high scatter within the individ- 
ual orders can be perceived. 




4500 4400 4300 4200 4100 4000 3900 3800 3700 3600 
observed wavelength [A] 

Fig. 7. All 42 lines with their radial velocity against the observed wave- 
length. Different orders are distinguished by colors and symbols. 



In the following, we consider a possible correlation of the 
line position uncertainty and its relative position within an or- 
der. Spectral distortions within the spectral orders have been in- 
vestigated at the Keck/HIRES spectrograph by comparing the 
ThAr wavelength scale with a second one established from 12- 
cell observations of a bright quasar by Griest et al. (1201 Oi l. In 
the wavelength range ~ 5000 - 6200 A covered by the iodine 
cell absorption they found both absolute offsets which can be as 
large as 500 - 1000 m s"' and an additional saw-tooth distortion 
pattern with an amplitude of about 300 ms"'. The distortions 
are such that transitions at the order edges appear at different 
velocities with respect to transitions at the order centers when 
calibrated with a ThAr exposure. This would introduce relative 
velocity shifts between different absorption features up to a mag- 
nitude the analysis with regard to A/z/yU is sensitive to. 



Whitmore et al. (12010b recently repeated the same test for 
UVES with similar finding though the saw-tooth distortions 
show slightly reduced peak-to-peak velocity variations of ~ 200 
m s '. The physical explanation for those distortions is not yet 
known, so it is still to be examined whether the deviations are the 
same at other wavelengths or depend on the specific exposure. 

The available solar atlas can be used to check UVES in- 
terorder distorsions as suggested in Molaro & Centurion ( 1201 11 1. 
For this purpose UVES observations of the solar spectrum re- 
flected by the asteroid Iris were taken on Sep 2009 with a re- 
solving power R s; = 85000. 

The differences of the positions in the UVES spectrum and 
the absolute positions of the solar atlas for 238 solar photo- 
spheric lines in the region between 500-530 run (Orders 121 to 
116) are shown in Fig [8] The schematic saw-tooth pattern de- 
tected by Whitmore et al. (2010) is also sketched in the figure. 
The stochastic distribution of the data points does not allow for 
any conclusion to be drawn about an underlying pattern. The 
observed dispersion is of 82 ms"'. Since the typical error in 
the measurement of lines in the UVES Iris spectrum is of =s 30 
m s^'and wavelength calibration residuals are of 25 m s"', there 
is an excess in the observed dispersion which suggests the pres- 
ence of local deviations in the UVES spectrum. The saw-tooth 
pattern detected by Whitmore et al. (2010) is not revealed by our 
test. 
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Fig. 8. Radial velocity shifts of individual lines Molaro et al. (201 1). 



Figure l9] sorts the observed lines according to their relative 
position within their order The origin of the abscissa reflects 
the central position within an order All observed orders are 
stretched to an identical scale and over-plotted for this purpose. 
Lines near -1/2 on the X-axis are positioned near the left rim of 
the order, and so on. The distribution of obtained radial velocities 
seems to show a certain periodic pattern. The blue curve shows 
a fitted cosine with an amplitude of 151 m s '. Considering the 
errorbars of the individual lines, this is no more than a slight indi- 
cation which supports the presence of local distortions resulting 
from a non perfect calibration. Systematic errors at the level of 
few hundred m s ' have been revealed also in the UVES data by 
comparison of relative shifts of lines with comparable response 
to changes of fundamental constants (Centurion et al. 12009 l l and 
Molaro et al (2012). Molaro et al (2012) suggest that these dis- 
tortions may originate from the block stiching of the CCDs. A 
CCD device is built-up by means of several sub-unit blocks with 
typical sizes of 512 pixels. The stitching of the blocks process 
produces misalignments of the order of few 0.01 of the pixel 
size in the block conjunctions. The ThAr has not enough lines to 
follow these imperfections which are therefore flattened in the 
pixel-to-wavelength conversion by a low order polynomial re- 
sulting into the observed spectral distortions. 
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Fig. 9. All 42 lines with their radial velocity against their relative po- 
sition within their order. A cosine fit with an amplitude of 151 m s^'is 
shown in blue. 



7. Discussion 

The result of Ayu/yu = (4.3 ± 7.2) x 10"^ we obtained is consistent 
with no variation of /i between Zabs = 3.025 and z = 0. The null- 
resuh is in agreement with recent publications on A/i/yu by King 
et al. ( [20m ) and Weerdenburg et al. ( 120111 ) at z^bs =2.811 and 
Zabs - 2.059, respectively. However, the present work utilizes the 
line-by-line fitting method (as i.e. in Ivanchik et al. l2005l l in con- 
trast to the other works which applied a comprehensive fitting 
method (CFM). The H2 system in the spectrum of QSO 0347- 
383 has the particular advantage of comprising a single velocity 
component, which renders observed transitions independent of 
each other. For absorption systems with two or more closely and 
not properly resolved velocity components many systematic er- 
rors may influence distinct wavelength areas. The CFM fits all 
H2 components along with additional H/ lines and handles an 
artificially applied A/u/fi as free parameter in the fit. The best 
matching A/z/yU is then derived via the resulting x^ curve. The 
CFM aims to achieve the lowest possible xl via additional ve- 
locity components. In this approach the information of individ- 
ual transitions is lost though since merely the overall quality of 
the comprehensive model is judged. 

In Weerdenburg et al. (2011 ) the number of velocity com- 
ponents is increased as long as the composite residuals of sev- 
eral selected absorption lines differ from flat noise. The residuals 
therein do not take into account the known inaccuracy of the es- 
timated flux error (see, i.e., Wendt & Molaro l20lll King et al. 
I20m i. 

As pointed out by King et al. ( 1201 11 1, for multi-component 
structures with overlapping velocity components the errors in 
the line centroids are heavily correlated and a simple x^ regres- 
sion is no longer valid. The same principle applies for co-added 
spectra with relative velocity shifts. The required rebinning of 
the contributing data sets implements further autocorrelation of 
the individual 'pixels'. 

The uncertainties of the oscillator strengths ft that are stated 
to be up to 50% (Weerdenburg et al. 120111 1 might further af- 
fect the criteria for additional velocity components. The method 
of CFM was applied for QSO 0347-383 by King et al. ([20081 . 
Section |4] and in particular Figure [3] reveal some of the men- 
tioned difficulties. The approach to fit individual lines with com- 
mon physical properties which was applied here allows us to 
carry out an error analysis which reproduces the impact of dif- 



ferential shifts within the spectral orders. This yields a higher 
transparency of the error-budget for individual lines but at the 
possible cost of larger scatter resulting in a slightly larger error- 
estimate. For the small number of H2 lines observed in the spec- 
trum of QSO 0347-383 we prefer the method applied in this 
work. The immediate advantage is that we are not forced to 
estimate the different systematic errors based on assumptions. 
Instead the true limiting error can be gathered directly from the 
data distribution and we are further able to attribute it to different 
sources. 

The new set of UVES observations of QSO 0347-383 this 
analysis is based has been taken with special care aimed to im- 
prove the measurement of Afif/u in a robust manner. In partic- 
ular the observations have been taken with higher resolution, a 
1x1 binning and calibration lamp spectra in direct combina- 
tion with the main exposures. These boost the precision of the 
analysis roughly by a factor two with comparison with Wendt 
et al. (12008b . We have shown that at the current level, calibra- 
tion issues become the dominant source of error. In addition to 
positioning errors, which are related to the signal-to-noise-ratio 
of the data, we observe for the first time inter order distortions 
which seem to be of the same order of magnitude of the uncer- 
tainties in the line positions. The conclusion is that we do not 
detect change in the value of // to 1 part in 10^ over a time span 
of 11. 5 Gyr, which is approximately 80% of the age of the uni- 
verse. 

High resolution data and attached calibration spectra are the 
key to understanding and handling the systematics which limit 
the precision of Afif/j measurements. It is important to fully con- 
trol the analysis of individual absorption systems and to mini- 
mize the errors involved wherever possible before extending the 
Ayu/ju analysis to multiple systems. Different characteristics of 
individual absorbers tend to get lost while not all errors of the 
measurements are likely to average out. 

Until new ways of wavelength calibration such as optical 
laser frequency combs (see i.e. Steinmetz et al. I2008I I are in- 
stalled for large optics, the data at hand is of the best quality 
available. High resolution VLT-data with special care regarding 
the calibration frames allows for the best precision that can be 
reached nowadays. In the context of Afi/fi measurements, opti- 
cal spectra of H2 at high redshifts are still of high importance to 
complement high precision determinations of Aju/// in the local 
universe via observations in the radio regime. 

Acknowledgements. We are thankful for helpful discussions on this topic with 
S.A. Levshakov and D. Reimers. 
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Table 3. Line parameters for the observed three rotational levels of H2 . 



J b[kms-'] logN 

~l 1.34 ±0.04 14.36 ± 0.02 

2 1.39 ±0.32 13.72 ±0.02 

3 2.02 ±0.05 13.91 ±0.01 



Table 4. List of 42 H2 lines fitted in QSO 0347-383. Radial velocities given relative to a redshift of z = 3.02489817. 



Line ID 


Lab. wavelength [A] 


Obs. wavelength [A] 


Pos. error [A] 


Velocity [km s ' ] 


Vel. error [km s '] 


K factor 


LlPl 


1094.0519 


4403.4530 


0.0026 


0.364 


0.174 


-0.003 


LlRl 


1092.7324 


4398.1360 


0.0015 


-0.039 


0.105 


-0.001 


L2R1 


1077.6989 


4337.6252 


0.0016 


-0.203 


0.113 


0.005 


L3P1 


1064.6053 


4284.9315 


0.0012 


0.247 


0.083 


0.010 


L3R1 


1063.4601 


4280.3151 


0.0011 


-0.242 


0.074 


0.011 


L4P1 


1051.0325 


4230.2996 


0.0021 


0.075 


0.148 


0.016 


L4R1 


1049.9597 


4225.9832 


0.0016 


0.165 


0.111 


0.016 


L5P1 


1038.1570 


4178.4767 


0.0018 


0.024 


0.126 


0.021 


L5R1 


1037.1498 


4174.4220 


0.0016 


-0.031 


0.113 


0.021 


L7R1 


1013.4369 


4078.9777 


0.0014 


-0.200 


0.104 


0.030 


L8R1 


1002.4520 


4034.7650 


0.0015 


-0.168 


0.109 


0.034 


L9P1 


992.8096 


3995.9617 


0.0017 


0.308 


0.131 


0.037 


L9R1 


992.0163 


3992.7595 


0.0015 


-0.389 


0.116 


0.038 


LlOPl 


982.8353 


3955.8147 


0.0019 


0.204 


0.146 


0.040 


LlORl 


982.0742 


3952.7519 


0.0017 


0.236 


0.127 


0.041 


L13P1 


955.7082 


3846.6281 


0.0033 


0.004 


0.255 


0.048 


L13R1 


955.0658 


3844.0424 


0.0024 


0.002 


0.187 


0.048 


L14R1 


946.9804 


3811.5054 


0.0025 


0.454 


0.198 


0.046 


WORl 


1008.4982 


4059.1012 


0.0017 


-0.093 


0.125 


-0.005 


WlQl 


986.7980 


3971.7622 


0.0016 


0.043 


0.118 


0.005 


W2Q1 


966.0961 


3888.4364 


0.0014 


-0.165 


0.107 


0.014 


L4P2 


1053.2842 


4239.3646 


0.0034 


0.208 


0.237 


0.013 


L4R2 


1051.4985 


4232.1695 


0.0024 


-0.350 


0.173 


0.015 


L5R2 


1038.6902 


4180.6199 


0.0032 


-0.162 


0.228 


0.020 


W0Q2 


1010.9385 


4068.9194 


0.0022 


-0.383 


0.160 


-0.007 


W0R2 


1009.0250 


4061.2231 


0.0038 


0.027 


0.278 


-0.005 



L2P3 1084.5603 4365.2445 0.0033 -0.008 0.227 -0.001 

L2R3 1081.7113 4353.7758 0.0027 -0.135 0.183 0.001 

L3P3 1070.1408 4307.2077 0.0025 -0.011 0.176 0.005 

L3R3 1067.4786 4296.4891 0.0022 -0.247 0.157 0.007 

L4P3 1056.4714 4252.1847 0.0027 -0.350 0.188 0.011 

L4R3 1053.9761 4242.1506 0.0018 0.307 0.125 0.013 

L5R3 1041.1588 4190.5564 0.0024 -0.137 0.168 0.018 

L6P3 1031.1927 4150.4420 0.0031 -0.253 0.226 0.020 

L6R3 1028.9866 4141.5628 0.0023 -0.256 0.163 0.022 

L8P3 1008.3861 4058.6547 0.0027 0.252 0.196 0.028 

L10R3 985.9628 3968.4022 0.0038 0.183 0.287 0.036 

L12R3 967.6770 3894.7970 0.0023 -0.332 0.180 0.037 

W0Q3 1012.6796 4075.9375 0.0021 0.381 0.157 -0.009 

W1R3 987.4487 3974.3837 0.0033 0.251 0.250 0.005 

W2Q3 969.0493 3900.3246 0.0022 -0.019 0.168 0.011 

W3P3 951.6718 3830.3853 0.0037 0.236 0.289 0.021 



